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It is known that the viscosity of a dilute gas can be derived by using kinetic 
theory. We present here a new derivation by using two entropy production principles: 
the steepest entropy ascent (SEA) principle and the maximum entropy production 
(MEP) principle. The known result is reproduced in a similar form. 
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INTRODUCTION 

Viscosity (or shear viscosity), 77, is a coefficient that appears in Newton's law of viscosity, 



r zx = V— , (1) 



where t zx is the shear stress and ^f- is the velocity gradient. For a dilute gas, the viscosity 

can be derived by using kinetic theory [11], with the result 
T3 ! 
C ■ 1 . 

f] = ~P vX i ( 2 ) 

where p is the density, v is the mean speed and A is the mean free path. Using kinetic 
theory to derive the viscosity is simple. We know that, for a given velocity gradient ^-, 
the x-component of the molecular momentum increases at rate M^- in the z direction; 
here M is the mass of a molecule. We also know that in a given direction, on average a 
molecule needs to travel X z to make a collision. So one may study a layer of a dilute gas with 
(r^ I thickness 2X Z in the z direction and view the dilute gas as two layers, each having height 

X z . Two such layers of dilute gas will exchange molecules and momentum. For each pair of 
molecules to be exchanged, the amount of x- momentum to be exchanged is twice \ Z M^-. 
During unit time, across unit cross-sectional area, the number of molecules to be exchanged 
is nv z , where n is the number density. The corresponding x- momentum to be exchanged 
is nv z X z M^-, which is actually the shear stress t zx . By using Eq. (JI|) and the relations 
\ z = 4=A, v z = -y=v and p = nM, the viscosity Eq. (j2]) is obtained. 

We shall present in this letter a new derivation of the viscosity of a dilute gas, by analogy 
with a derivation of the thermal conductivity of a dilute gas [2] . The new derivation will 
use the entropy production. Here, two aspects of the entropy production will be used: 
(1) for a given state, it evolves in the direction of the steepest entropy ascent; (2) among 
many candidates, the actual steady state has the maximum entropy production. We call 
the first the steepest entropy ascent (SEA) principle and the second the maximum entropy 
)roduction (MEP) principle. The name of "SEA" was originally introduced by Beretta [3- 
to postulate that a system always evolves to the direction having the steepest entropy 
ascent. The MEP principle has been used by Paltridge [8l-ll0l] to study the Earth's climate by 
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postulating that the steady state of the atmosphere has the maximum entropy production. 



Though in their original proposed form, the SEA principle and the MEP principle overlap 
and one even can be alternative to the other [7|, we use here their names to emphasize two 
separate aspects of the entropy production. There exist many studies of entropy production 



theory itself; see, for example |lll421j. This paper provides an application. 

In order to derive the viscosity of a dilute gas, we shall study a dilute gas confined between 
two parallel plates where the bottom plate is stationary and the top plate is subject to a 
constant shear stress. Let the dilute gas be subject to three restrictions: (1) the density p 
or n is constant, and so we can write the number of molecules as iV = nAAz, where A is 
the cross-sectional area and Az is the height of the dilute gas; (2) N is large, N ^> 1, and 
so we can apply some approximations; (3) the shear stress is small, and so only laminar flow 
appears, and the velocity gradient is small. Such a dilute gas will carry a steady velocity 
gradient from whose expression the viscosity can be extracted. The expression of the steady 
velocity gradient will be determined by using the two entropy production principles. To do 
that, we shall first obtain the entropy of a dilute gas with respect to the velocity gradient. 



ENTROPY OF A DILUTE GAS HAVING HEIGHT A 2 



In this section, we study a dilute gas having height X z and estimate its entropy with 
respect to the velocity gradient. As Fig. [TJand Fig. [2] show, for each molecule, we only need 
to study how it moves in the x direction, given that the velocity gradient is ^f. So we may 
simply consider that each molecule moves either in the +x direction or the —x direction at 
the same speed v z . 




FIG. 1: The molecular speed distribution of a dilute gas in the equilibrium state. The left graph 
shows only the molecular speed in the x and z direction. The middle graph shows only the 
molecular speed in the x direction. The right graph shows only whether a molecule moves along 
the +x direction or the —x direction. The right graph can also be viewed as that a molecule 
moves either in the +x direction or the — x direction with the same speed v z , which satisfies 

— — — V 
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Let us study a macroscopic state in which k molecules out of N molecules move in the 
+x direction. The number of corresponding microscopic states is Q(k) = C^, and by using 
Boltzmann's entropy formula the corresponding entropy is 

iV' 
S(k) = k B In fi(fc) = k B In C k N = k B In - (3) 



k\(N-k)\ 



FIG. 2: The molecular speed distribution of a dilute gas carrying a velocity gradient -^-f . The 
dilute gas is on a stationary plate, its height is A 2 , and it consists of N molecules of which k 
molecules move in the +x direction. It can be viewed as a fluctuation of the equilibrium state 
in Fig. Q] for which k « y, but we have here k = y + -4^-37 • For a given k, the number of 
corresponding microscopic states is C 1 ^ and the entropy is S(k) = /cslnC^-. 



Given the velocity gradient yf- being small, k is comparable to y, k ~ y. Given iV being 
large, Stirling's approximation can be used, In iV! m NlnN — N and In A;! w klnk — k. Then 



2 fc 

by using Taylor's series approximation ln(l + x) ~ x — y- , where x — - — 1 



2-, we get 

2 
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(4) 



where So is the entropy of the equilibrium state. Given that k molecules move in the +x 
direction, the number of molecules that move in the — x direction is N — k. Then the 
macroscopic x-velocity of the whole dilute gas is 

(2k - N)v z 



N 



(5) 



The same velocity can also be obtained as the average macroscopic directed flow velocity 

dv x X z 



dz 2 ' 
The two expressions are about the same quantity and so they are equal, which leads to 

N NX Z dv x 

2 Av z dz ' 
Then by using N = nAX z , we may rewrite Eq. (jl]) as 



(6) 



(7) 
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ksnAXl ( dv. 



dz 8v 2 \ dz 

Actually, the entropy is still lower, as discussed in the next section, and we write 

^,dv x , „ kunAXl ( dv x 

^(-7 £ ) ~So- -^—^ -f- 

dz Ivi \ dz 



(8) 



(9) 



This is the entropy for a dilute gas that has height of only A 2 



DISCUSSION 

In our study, we have chosen the height of the dilute gas as X z . By choosing height A^, 
a molecule can collide both with other molecules and the bottom plate, which is stationary. 
When a molecule collides with the bottom plate, its average velocity remains about zero; 
when molecules collide with each other, a velocity gradient may appear as a result of a 
velocity distribution fluctuation. 

We have written the entropy of such a dilute gas as Eq. Q, which is Eq. (jSJ) plus two 
corrections. 

(1) Eq. (JSJ) is derived from Eq. (J3J). But given a k, the corresponding velocity gradient 
does not necessarily appear, it only appears with an additional feature. The velocity gradient 
is distributed uniformly across the whole dilute gas, although local fluctuations also exist. 
Thus the dilute gas can be viewed as two sub-layers where each sub-layer can be analyzed 
in the same way as the original layer, see Fig. [3J By recursively going through this process 
and neglecting the local fluctuations, Eq. (JHj) becomes 
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FIG. 3: Additional negative entropy arising due to the uniform velocity gradient distribution. 



The left graph shows a dilute gas having height X z and carrying a velocity gradient 



dv x 
dz 



The 



middle graph is for the same dilute gas but has been boosted to the rest frame. The right graph is 
also for the same dilute gas in the rest frame but has been viewed as two dilute gases each having 
height A z /2. Then the discussion about Fig. [2] applies to each of them, except that the height now 
becomes X z /2. Thus, additional negative entropy arises and Eq. (|5J) becomes Eq. (fTU]) . 

(2) Given a velocity gradient =¥*■, the molecular x— speed distribution is not the same 
as the equilibrium state, even in the center-of-mass frame. Think of a dilute gas that can 
be viewed as many layers each having height X z ; any two adjacent layers will have different 
macroscopic velocities, see Fig. H] and the corresponding kinetic energy is 



2 x -pAX 2 



dv x X z 
~dz^2 



pAX 3 z ( ' dv x 

dz 



(11) 



This energy, when the velocity gradient relaxes, will become heat and lead to an additional 
entropy 



where T is the temperature. By using p = nM and v = \/^fj-, we write this entropy as 



2kBnAX 3 ( dv x \ 2k B nAX 3 ( dv x ^ 
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(13) 



By this amount, the entropy of the dilute gas is further reduced and so Eq. (ITU]) approxi- 
mately becomes Eq. (jUJ). 




dv x A 2 
dz 2 


• 








< — < 

• 



d^ 2 



FIG. 4: Additional negative entropy of a dilute gas arising due to the large height. The left graph 
shows a dilute gas carrying a velocity gradient ^~, its height is 2A Z . The middle graph is for the 
same dilute gas but has been boosted to its rest frame. The right graph is for the same dilute gas, 
but has been viewed as two dilute gases, each of which has height X z and is in its rest frame, with 
the two dilute gases moving in opposite directions with velocity -jjf-f- For each dilute gas, this 
velocity leads to a kinetic energy p g * (-gj-) and an additional negative entropy P 8T Z (-^f ) • 



ENTROPY PRODUCTION AND VISCOSITY OF A DILUTE GAS 



In this section, we start with the calculation of the entropy production of a dilute gas. 
The dilute gas is confined between two plates, with the bottom plate stationary and the 



top plate subject to a constant shear stress t zx . Across the dilute gas, a steady velocity 
gradient will be established. Its expression shall be determined by using the SEA principle 
in conjunction with the MEP principle. To do that, for a given velocity gradient, we need to 
examine two specific possibilities of how it evolves and calculate the corresponding entropy 
production. 

One possibility is that the given velocity gradient is steady. For this possibility, the heat 
generated during unit time interval is 

r«AAz^, (14) 

where A is the cross-sectional area of the dilute gas and Az is the height. This heat will be 
absorbed by the dilute gas. The entropy of the dilute gas (or the environment if the heat is 
transfered to the environment) will increase. The corresponding entropy production is 

"=^t£- (15) 

where T is the temperature of the dilute gas. 

The other possibility is that the given velocity gradient ^ starts to relax at the maximum 



rate, 



dz 

£„(-!). (16 > 



where r is the minimum possible relaxation time, r is approximately X/v. To see why 
the relaxation is exponential and why r « X/v, let us study an example of a dilute gas 
that has height 2X Z and is confined between two plates moving in opposite directions at 
the same speed. A velocity gradient will be established across the dilute gas. Suppose that 
the two moving plates are suddenly detached from the dilute gas, the velocity gradient will 
relax gradually at the maximum rate. Let us view the dilute gas as two layers each having 
height X z , and study their relative velocity. The relative velocity will relax as a result of the 
momentum exchange between the two layers, but the momentum exchange also depends on 
the relative velocity and is proportional to it. Thus, the smaller the relative velocity is, the 
smaller the momentum exchange rate becomes; the smaller the momentum exchange rate 
becomes, the slower the rate that the relative velocity will relax. As a result, they both 
decrease exponentially. Thus the velocity gradient relaxes exponentially too. Concerning 
r, a molecule needs to take time X/v to transfer a momentum from one layer to another, 
and the mean value of the momentum that it transfers is comparable to the mean molecular 
momentum difference between the two layers. Thus, within a few transfers for each molecule, 
the velocity gradient will relax. So the relaxation time r is approximately X/v, but a little 
bit larger. What we have studied here is a dilute gas that has height 2X Z . For a dilute gas 
that is thicker, it can be viewed as many layers that each has height X z . Then, between any 
two adjacent layers, the discussion is the same. But each layer except for the two end layers 
has two adjacent layers and thus all layers are correlated. This makes the relaxation time 
larger. But we are only concerned with the instant at which the velocity gradient starts to 
relax; at that instant, the relaxation correlation has not yet taken effect, and thus, in the 
limit t — > 0, r should be the same for a dilute gas having any height. 
Given the velocity gradient Eq. ( TT5|) . Eq. (J5J) leads to 

m = So _!^^vr _ 2 iy (17) 



and the corresponding entropy production is 



dt 
which, in the limit t — > 0, leads to 



aUi j_m__j_^i(^) exp ( 2 1, 
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This is the entropy production of a dilute gas that has height A^. For a dilute gas thicker, 
in the limit t — >• 0, it can be viewed as many such layers and the entropy production is then 

Az ksnAXf ( dv x \ . . ksnX 2 ( dv x \ , . 

'= K^fiii) =AAz —(n) ■ (20) 

This is because at the instant when the velocity gradient starts to relax, any two non- 
adjacent layers should have no relaxation correlation established yet. The reason why we 
use the limit t — Y is that we only need to know whether or not the given velocity gradient 
would start to relax. The relaxation process itself is irrelevant. 

The entropy productions of the two possibilities are plotted in Fig. They are equal 
when ^ = (^f) <• We shall use first the SEA principle and then the MEP principle to 
determine that the steady velocity gradient must be (^-) A - According to the SEA principle, 
for a velocity gradient ^j > (^j) A , it cannot be steady because the entropy production 
would not be the maximum. (But it would not start to relax at the maximum rate either; 
actually, the velocity gradient will start to relax at a smaller rate until ^ = ( ^ ) *, and 
during the entire relaxation process the shear stress will continue to generate heat and 
entropy.) Thus the steady velocity gradient must be in the range < ^ < (^) A - Then 
according to the MEP principle, we see that the steady velocity gradient must be (^fj-K 
because, among all candidates, it is the one that is associated with the maximum entropy 
production. 

Another way to think about Fig. [5] is that the two entropy productions compete with 
each other. Given a dilute gas filled between a stationary bottom plate and a moving upper 
plate, the dilute gas itself tends to relax toward the equilibrium state, while the upper 
moving plate tends to move faster and faster given the shear stress being constant. But 
the two tendencies are in opposite directions and compete with each other, until at last one 
balances the other and the velocity gradient becomes steady. Even so, the relaxing tendency 
still presents. And its strength can be measured by the entropy production shown in Eq. 
(120]) . So, for a given velocity gradient to be steady, it must at least overcome that relaxing 
tendency by having an equal or larger entropy production. Thus, its upper limit is obtained 
from Fig. 0as (j) A - At the same time, the upper plate still tends to move faster, and its 
strength should be measured by another entropy production, which must be overcome too 
if a velocity gradient is to be steady. Thus the lower limit of the steady velocity gradient 
should also exist. And perhaps the upper limit equals the lower limit and they are the 
steady velocity gradient. If so, the steady velocity gradient can be obtained on the basis of 
the SEA principle alone. 

By calculating (^f) A in Fig. [5j the steady velocity gradient is obtained as 

dv x tv 2 

dz k B TnX 2Tzx - [ } 
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FIG. 5: Entropy production with respect to small velocity gradient. The solid curve is the entropy 
production when the velocity gradient is steady. The dashed curve is the entropy production when 
the velocity gradient starts to relax at the maximum rate. {-&[-) A, the velocity gradient at point 
A, is the determined steady velocity gradient. For a velocity gradient -^ > (-j^-)a, it cannot be 
steady, because in this range the solid curve is not the highest, at least it is lower than the dashed 
curve. The steepest entropy ascent (SEA) principle indicates that a state can occur only if it is 
associated with the maximum entropy production. Thus the steady velocity gradient must be in 
the range < -^f < (-^f )a- Subsequently, according to the maximum entropy production (MEP) 
principle, we know that the steady velocity gradient must be {-§§-) A, because among all candidates 
it is associated with the maximum entropy production. 



So we have 



k B Tn\ 2 dv x 



1 ZX _o , 

tv z dz 
Comparing it with Eq. (0Q), the viscosity of the dilute gas is extracted as 

k B Tn\ 2 



V 



TV 2 



By using relations 



'8k B T %r A 

r = \l — — , p = nM, r ~ -, 

7TM V 



we may write it as 

which is comparable to the kinetic theory result 



i] = —pvX, 



(22) 

(23) 

(24) 
(25) 



CONCLUSION 

We have derived the viscosity of a dilute gas by using a new approach. The new approach 
is based on the fact that there exists an entropy production competition between at least 
two possibilities of how a given velocity gradient evolves. One possibility is that the velocity 
gradient is steady. The other possibility is that the velocity gradient relaxes exponentially 
at the maximum rate. Then, by applying the steepest entropy ascent (SEA) principle, we 
have obtained the upper limit of the steady velocity gradient. Subsequently, by applying 
the maximum entropy production (MEP) principle, we have identified the actual steady 
velocity gradient from whose expression the viscosity is extracted. The viscosity extracted 
is comparable to the result of kinetic theory. 
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